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' We develop a thorough description of neutrino oscillations in a magnetized protoneutron star, 
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J>^' based on a resonance layer for neutrinos with different momentum directions. We apply our 
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approach to the calculation of the asymmetry in the neutrino emission during the birth of a 

0^ ' neutron star and the pulsar acceleration in the case of an active-sterile neutrino resonant conversion. 
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protoneutron star, for sterile neutrino masses of the order of KeV and small mixing angles. 
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I. INTRODUCTION 



The peculiar motion of pulsars is one of the challenging problems in modern astrophysics. 
Observations show that young pulsars move with velocities of 200 — 500 km s^^ on average, 
much larger than the mean velocity of ordinary stars in our Galaxy (~ 30 km s-^)^. This 
is strong evidence that pulsars receive a kick during their birth in supernovae explosions, 
but the physical origin of such an impulse remains unclear. A natal kick could be a result 
of global hydrodynamical perturbations in the pre-coUapse core of the progenitor star or 
an anisotropic neutrino emission induced by the strong magnetic fields in the protoneutron 
starj^. A significant fraction of the pulsar population has velocities over 1000 km 
which are difficult to account for assuming hydrodynamical kick mechanisms, as numerical 
simulations of supernova explosions indicate. 

More than 99% of the gravitational energy released in the collapse is emitted in the 
form of neutrinos and antineutrinos. A small asymmetry (~ 1% ) in the momentum taken 
away by neutrinos could produce an acceleration consistent with the measured pulsar veloc- 



ities. A particle physics exr 



anation of the origin of this asymmetry has been proposed by 



Kusenko and Segre (KS)P]|^, based on adiabatic flavor (z/g < — > ^ij.,t) oscillations in matter, 
in the presence of a magnetic field. In a magnetized stellar plasma there are anisotropic 



contributions to the neutrino refraction index! 



, and thus the resonant transformation 



of neutrinos moving in different directions respect to the magnetic field occurs at different 
depths within the protostar, producing a momentum flux asymmetry. With the standard 
(active) neutrinos, the mechanism works when the resonance region lies between the sur- 
face of the two neutrinospheres, but it requires an exceedingly high square mass difference 
(Am^ > 10^ eV'^), in conflict with the existing limits for standard neutrinos. 

A variant of the idea of KS can be implemented in terms of matter neutrino oscillations 
between active and non interacting (sterile) states Q], in which case the above mentioned 
limitation can be avoided. More recently it has been shown that an asymmetric off-resonant 
emission of sterile neutrinos is also a feasible mechanism^]. In both cases, the range of 
required oscillation parameters (mass and mixing angle) overlap with the allowed region for 
sterile neutrinos considered as the cosmological dark matter. Alternative neutrino oscillation 
mechanisms for the pulsar kick have been proposed in relation to the possible existence of 
transition magnetic moments j9i], non orthonormality of the flavor neutrinos jlo|. or violations 
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of the equivalence principle The asymmetric neutrino emission could also produce a 
short gamma-ray burst |l2l| and gravity waves, whose detection can help to understand pulsar 
kicks |l3l]. 

In a previous article [ij] we reexamined the viability of the KS mechanism for resonant 
active-active neutrino oscillations, and emphasized the relevance of the geometrical 
dependence of the energy flux to produce an anisotropic neutrino emission. Following previ- 
ous works on the subject, we use the concept of an induced deformation by the magnetic field 
of an averaged resonance surface, which acts as an effective u^^^- neutrinosphere. In this work 
we analyze the active-sterile transformations incorporating a more realistic description in 
terms of a spherical layer where the resonant conversion of neutrinos with momentum k mov- 
ing in all directions takes place. A geometrical expression for the neutrino flux anisotropy 
is derived in terms of the ratio of the thickness of the resonance shell, which depends on the 
intensity of the magnetic field and its average position. To quantify the effect we apply our 
results to two simple models for the protostar atmosphere. 

The paper is organized as follows. In Section II we analyze the resonant conversion of 
active neutrinos into sterile neutrinos within the protostar atmosphere in the presence of a 
magnetic field, and we derive an expression for the thickness of the resonance region as a 
function of the magnetic field and the medium properties. In Section III, on the basis of 
the reasonable assumption that the neutrino flux is in the diffusive regime, we evaluate the 
asymmetry of the momentum flux taken out by neutrinos. In Section IV concrete results 
are obtained in the context of two specific models. Section V contains some comments and 
concluding remarks. 

II. RESONANT ACTIVE-STERILE TRANSFORMATION IN A MAGNETIZED 
MEDIUM 

Within the neutrinosphere neutrinos are trapped due to the opacity of the medium. This 
region can be considered as limited by a spherical surface with a radius R,, , , defined by the 
condition that above the surface neutrinos suffer less than one collision Il5l 




(1) 
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where X^^ is the mean free path of the neutrino of flavor I = e, ^, r. We consider the proto- 
star atmosphere as constituted by a degenerate gas of relativistic electrons and a classical 
nonrelativistic gas of nucleons. The main contributions to the neutrino opacity in such a 
medium come from the neutral-current scattering reactions jQ 



i^i + n ^ Ui + n, 
ui+p ^ 1^1 +p, 

and the charged- current absorption reaction 

Uf. + n ^ e~ + p. 

The cross sections for the reactions Q and Q are 



8 sin^ Ow 



'I - 2sin2 Ow) 



(2) 
(3) 

(4) 

(5) 
(6) 



1 + 3^/1 

where k is the neutrino momentum, sin^ 6w ~ 0.23, and qa — 1-26 is the renormalization of 
the axial-vector current of the nucleons. For the absorption reaction we have 



Q 

CTabs = 4(T„ ( 1 + — 



Q 



1/2 



(7) 



with Q = uin — mp denoting the difference between the neutron and the proton masses. 

In terms of thermally averaged cross sections, the mean free paths for the electron and 
muon (tau) neutrino are given by 



X^^ = Nn {{(Tabs) + {cTn)) + Np{ap) = K^^^k^p, 
K^r = ^n{(Tn) + Np{(Tp) = Ku^^^Pp, 



(8) 
(9) 



where N^^p are the neutron and proton number densities, p is the baryonic density, and 
are the neutrino opacities. In a protoneutron star Nn — 0.9 (A^„ + Np) and k ^ Q, which 
yields 



and 



= 3x 10"^^ MeV-^, 
n = 0.7 X 10-2^ MeV-^. 



(10) 

(11) 
(12) 



Eq. (fTUj) implies that R^^ > -Ry^_^, as can immediately be seen from Eq. ((T)). 

In the presence of a magnetic field, the effective potential for an active neutrino propa- 
gating in a plasma of electrons, nucleons, and neutrinos can be written asjsl [l7|: 



K, = K, - ec^, k ■ B 



(13) 



Ne+[^ -2sm^ 9w] (iVp - N,) - ^iV„ 



with the coefficients b,y^ and c^^ given by 

= V2Gf 
K = V2Gf 
c,, = 2V2Gf 

= 2V2Gp U fl + 2mp^) Gp - ^^2m„^C„ + C, 



+ K , 



( 1 + 2mp^j Gp - gA2mn^Gn - G^ 



(14) 

(15) 
(16) 

, . - (17) 

e / e J 

where e > is the proton charge, the anomalous part of the nucleon magnetic 

moments 

/s:„ = -1.91-^, Kp = 1.79-^, (18) 

and 



V2Gi 



a. 



d^p 1 d 



fe,n,p ■ 



(19) 



(20) 



''"'^ 2 J {2nf2EdE" 

In these equations N^. and Nu^{Np^) are the number densities of electrons and neutrinos 
(antineutrinos) of flavor /, respectively and fe,n,p are the distribution functions of electrons, 
protons, and neutrons. Since sterile neutrinos do not interact with the medium 

= v,^ = 0, 

which means that, unlike the oscillations among active neutrinos, the neutral current con- 
tributions to V^^ (Vpi) play a relevant role in the MSW equation governing the active-sterile 
oscillations in matter. 

For antineutrinos the potentials change in sign, i.e. Vp, = — Vj^,, and in general the 
transition probability for ui < — > Ug will be different for neutrinos and antineutrinos. As a 
consequence, their relative amounts will change, driving the potential locally to zero and the 



mixing angle in matter close to the one in vacuum throughout the protoneutron star|l8j. In 
this work we will assume that this effect is not important because the resonance conversion 
takes place within a thin layer. Under such conditions the kick mechanism happens as 
described in Ref. jsj. 

For relativistic degenerate electrons Cg = — /ie/Svr^ and for classical nonrelativistic nucle- 

ons Cn,p = -Nn,p/8Tm 

n,p 

, which gives 



-4^20^ ( ^ + + tS ) . (22) 

... . we ... . He. . . .e W„ _ .a . - ,3.^., 
is the electron chemical potential. It should be noticed that in Ref. j2| Cj^^ vanishes because all 
the components of the stellar plasma are assumed to be degenerated, and the contributions 
due to the anomalous magnetic moment of the nucleons are not included. These expressions 
have been derive in the weak-field limit , i.e. B ^ /ig/2eQ|. More general features of 
the neutrino propagation in magnetized media , in corporating the effect of strong magnetic 
fields, have been considered by several authors jl7|. 

The neutrino fractions are much smaller than the electron fraction Yg = ^e/ {Nn + Np) 
~ 0.1 and, in the context of our analysis, the contributions to K,, coming from the neutrino- 
neutrino interactions (denoted by hy^) can be neglected. With this approximation in mind, 
and taking into account that by electrical neutrality = Np, the coefficients in the effective 
potentials become: 

K = -h,^-^il-mYe)^, (23) 



3iV„ 


Pe 


20Tm„ 


167r2 




Pe 


20rm„ 


167r2 



r = — r- ~ 



where ?7e = 3, = 1 and r^^,^ = 1, = -1. 

As is shown by Eq. ()13p. the electron fraction plays a very important role in the case of 
active-sterile oscillations. This quantity can be estimated as follows. From the equilibrium 
condition for the reaction e~ + p ^ n + u^, taking ~ we have the relation fig + fip c:^ /i„ 
for the chemical potentials. Hence, Nn = 6"'-'^"^=-'/"^ Np for a non-relativistic nucleon gas at 
a temperature T and 

6 



where we used the fact that Q fie for typical densities in a protoneutron star. The above 
is an imphcit equation for Ye, which has to be solved numerically once p and T are known 
as a function of r. 

Now we use the effective potentials with bp^ and Cp^ given by Eqs. (f2H|) and (f2^ to discuss 
the effect of a magnetic field on the active-sterile matter neutrino oscillations. To make our 
analysis simpler, let us consider the case of two states (z/;, z/^) that are a mixture of the mass 
eigenstates ui and z/2. For a momentum k, a ui will undergo a resonant conversion into in a 
Us if the following condition is satisfied: 



■cos 20 = V 



2k 



where 6 is the vacuum mixing angle and Am^ 



•^l lr=_R{k) 



(26) 

ml is the difference between the 
squares of the neutrino masses. Assuming Am^cos20 > 0, if S = the resonant condition 
for the electron neutrino (antineutrino) is satisfied when Ye > 1/3 {Ye < 1/3). In addition, 
since Ye < 1 only the z/^^t (but not the z^/^.t) can go through the resonance region. In what 
follows we assume that within the protoneutron star Ye < 1/3, so that the resonant condition 
is verified only for antineutrinos. 

In a linear approximation, we can write 



R{k) = Rr{k) + 5{k) k-B, 



(27) 



where Rr{k) is the radius of the resonance sphere in the absence of the magnetic field. For 
a certain value of the magnitude of the momentum k, Eq. (jTTjl determines a spherical shell 
limited by the spheres of radii Rr{k) ± 6{k)B, where the resonance condition is verified for 
antineutrinos moving in different directions with respect to B. Using the expression (fT!?|) for 
Vp, in Eq. we get 



where 



and we have defined ^ 



n 



1 Cy 



--Rr{k) 



explicitly 



^In^f for any function g{r). From Eqs. 

2 n + Ye 



and 



;i - r^iYe) h-^ - ifiYehyl V 5rm„ Ar^^p 



(28) 

(29) 
we find 

(30) 



r=Rr(k) 



By means of Eq. ()25|1 . hy^ can be expressed in terms of h^^ and hj,'^ as follows 



h- 



yI — 



(31) 



The efficiency of the flavor transformation depends on the adiabaticity of the process. 
The conversion will be adiabatic when the oscillation length at the resonance is smaller than 
the width of the enhancement region. This imposes the condition 



with bp^ given by Eq. ((221) • 

III. NEUTRINO MOMENTUM ASYMMETRY 

The energy flux in the atmosphere is radial and determined by the diffusive part of 
the neutrino distribution function. Consequently, in previous works on the subject the 
resonance condition has been evaluated in terms of an average neutrino momentum (k) = 

7 4 

i8(3C(3) '^ ~ 3.15T pointing in the radial direction. This procedure leads to the concept 
of a deformed average surface of resonance, which acts as the boundary of an effective 
neutrinosphere. In a more detailed description, at each point of the neutrinosphere we have 
neutrinos with momentum k pointing in all directions. Hence, the resonance condition as 
given by Eq. actually deflnes a spherical surface for each value of k ■ B, which acts 

as a source of sterile neutrinos moving in the k direction. Accordingly, as we mentioned 
above, for a given k the resonant transition takes place within a spherical shell instead of 
a deformed spherical surface. We will discuss the situation where the resonance shell lies 
totally within the i/;-neutrinosphere. If the mixing angle is small and the adiabatic condition 
of Eq. is satisfled, then the ui are almost completely converted into z/^. The outgoing 
sterile antineutrinos freely escape from the protostar, while the ingoing ones propagate 
through the medium until they cross again the resonance surface and are reconverted into 
ui, which are thermalized. In consequence, only those Ug going outward leave the star and 




(32) 



From Eq. ()13j) and neglecting the contribution of the magnetic fleld we get 




(33) 
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the resonance surfaces for each momentum direction act as effective emission semispheres. 
For k pointing in opposite directions the respective semispheres are at Rr{k) ± 6{k)kB and 
therefore have different areas, which generates a difference in the momentum carried away 
by the sterile neutrinos leaving the star in opposite directions. This difference adds up 
to building a nonvanishing asymmetry in the total momentum /C emitted by the cooling 
neutron star, which we now compute explicitly in terms of the whole neutrino distribution 
function. 

Within the neutrinosphere the ui satisfy the diffusion regime with a distribution 

function 

/. =r ^-T^k-Vf', (34) 

where P*^ = (l + e^^'^) ^ is the distribution function at equilibrium. We have again assumed 
that the neutrino chemical potential is negligible. The momentum flux emitted in the radial 
direction is 

dFr = k(fk. (35) 

Let us now adopt a reference frame (x, z) fixed to the star, where the magnetic field 
coincides with the 2;-axis, B = i?z. In addition, at each point within the spherical shell 
determined by Eq. (j^7j) we use a local frame (x', y', z') with f = z', such that k = cos O'i' + 
sin 9' cos 0'x' + sin 9' sin 0'y'- this way the total momentum carried away by the vi can 
be expressed as 

OO l-TT /•2tT 

k^dk sm9d9 sm9'd9' d^'R\k)f_k-r (36) 



i2nrJo Jo Jo " ^ 



^-l (37) 

The momentum emitted in direction of the magnetic field ks is calculated in a similar 
manner, with k ■ f replaced by k • B = cos 9 cos 9' — sin 9 sin 9' sin 0': 

kB = — dk k'Rr{k)6{k) r - 7, r^-^ , (38) 

where we have used Eq. (|27|) keeping at most terms that are linear in 6{k). 

In terms of the above quantities the z/; contribution to the fractional asymmetry in the 
total momentum becomes 



where the factor 1/6 comes from the fact that we are assuming that K. is equipartitioned 
among all the neutrino and antineutrino flavors. In order to evaluate the remaining integrals 
in Eqs. (I37p and (|38|) we need to know the explicit dependence on k of the functions Rr and 
6. Simple analytical results can be derived by replacing these functions by their values at 
the average momentum (k) . Proceeding in this manner we get 

48 y 5 3k,^p dr j ' ^ ' 

In these expressions the terms that depend on dT'^/dr come from the diffusive part of the 
distribution function. In the regime considered here they happen to be smaller than the 
first terms and we will discard them. This can be verified in the context of specific models 
like the one considered in the next section. Using this approximation, we finally obtain 

A/C 2_6_ 

which is four times bigger than the asymmetry we derived by assuming a single (deformed) 
resonance surface [l^ . This means that if Eq. were applied to the case of active-active 
oscillations, when the resonance shell is at the border of the electron neutrinosphere, we 
would obtain B c:^ 7.5 x 10^^ G for the intensity of the magnetic field required to produce 
the observed pulsar velocities. 

Up to this point we have considered the resonant oscillation between a sterile antineutrino 
and an active antineutrino of a given flavor. The observational evidence indicates that active 
neutrinos mix. Furthermore, it is reasonable to assume that the three active antineutrinos 
could be mixed with the sterile ones. Consequently, the situation can be more complicated 
than the simple pairwise neutrino oscillation already discussed. In general different mixing 
schemes will render gives different results for the magnetic field, but these results will differ 
in factors of the order of the unity. Therefore, to make an estimation of the magnetic field 
we will consider the simplest possible scheme, in which the mixing between active neutrinos 
is neglected and pairwise oscillation with the sterile ones are taken. In this context the total 
momentum asymmetry can be written 



IC 9 \ Rr- Rr- 
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In what follows we will make the calculation for the active-sterile conversion in the framework 
of this approximation by assuming that all the resonance layers are entirely within the i^^^,-- 
neutrinosphere. 



IV. NUMERICAL RESULTS 

Taking into account the total momentum emitted by the protostar, to have the required 
kick A/C//C must be of the order of 0.01. Inserting this value for the momentum asymmetry 
on the left hand side of Eq. (j43p we see that ~ 0.015/?^^,^ , which means that the thickness 
of the resonance shell is of the order of one hundred of meters. From Eqs. and 
using Rr- ~ — Rr, we find 

eB = 0.045 (44) 

Vp +2Vp ^ ^ 

The magnitude of the magnetic field needed to produce the kick can be evaluated from this 
result and Eq. ()30|1 . once the temperature and density profiles are known. In this section 
we will perform this calculation for two simple models of a protoneutron star atmosphere. 
One of them is the spherical Eddington modelQ] and the other is the one introduced in 
Ref. l21|, on the basis of a fitting of numerical results. 

n 

The Eddington modelj20|] gives a satisfactory description of a neutrinosphere locally ho- 
mogeneous and isotropic, with a conserved energy flux. In the case of a protoneutron star 
the results are in reasonable agreement with those obtained by numerical analysis. For a 
spherical geometry the model was developed in Ref. , where the reader is referred for de- 
tails. Here we summarize those features that are relevant for our analysis. The atmosphere 
of the protostar is described as an ideal nucleon gas in hydrostatic equilibrium 

P = ^, (45) 



dP GM(r 
= -p 

ar r 



(46) 



with M(r) = 4:71 Jq dr' r''^p{r'). The total neutrino momentum flux satisfies the diffusion 
equation 

nr) = -T^-^ (47) 
6b Kp[r) dr 
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where k = { k + 2k _ ^ ) = 3 x 10 MeV ^, and the flux conservation 

nr) = j^. (48) 

where is the luminosity of the protostar. Eq. is obtained from Eq. by integrating 
over the momentum. 

The above system of equations has no exact analytical solution. However, a good ap- 
proximate solution has been found in Ref. Q|. It allows us a consistent treatment, from 
which we obtain the following expressions: 

- - ^ • («) 

+ = (SO) 

with Ac = {9kLcPc) / i'^.TiT^ Rc) . Here pc and Tc are the density and the temperature at the 
radius of the core Rc, and a is a function of the temperature 

where = [GMcrrin) / {2XcTcRc) and Tg is the asymptotic value of the temperature. In Fig. 
1 we show the resulting profiles for p(r), T{r), and M(r), corresponding to a configuration 



like the one considered in Ref. [1^: M^. = Mq, Rc = 10 km, = 2 x 10^^ erg s~^, 
Pc = 10^'' g cm~^, Tc = 40 MeV, and Tg = 4.8 MeV. In the same figure, we plot Ye{r) as 
obtained by solving Eq. ()25|) at every point in the atmosphere. 
According to Eq. pO|) the electron neutrino mean free path is 

7) (t)' («^) 

The radius of the neutrinosphere is estimated by inserting the above expression into Eq. 
and evaluating the resulting integral numerically, the result is Rp^ ~ 2.6 Rc- This yields 
p{RpJ — 1.3 X 10^^ g cm~^, T{RpJ ~ Tg , and Ye{RpJ ~ 0.1 which are in good agreement 
with typical values for the electron neutrinosphere in a protoneutron star. In the same way, 
for the muon (tau) neutrino we get Rp^^^r — 2-2 -Re- 
using Eqs. PHj) and (07j), the magnitude of the ratio of the second and the first terms 
in Eq. pn|) is {LcKyi) / {^tRT^r"^), which varies between 10~^ and 10"^ for Rc ^ r ^ -^j^/^.r- 
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Therefore, the approximation done in writing Eq. ()42|) is justified and the intensity of the 
magnetic field required to produce the pulsar kick can be evaluated from Eq. (03}. As shown 
in Fig. 2, at the surface of the core — 8 x 10^^ G and B remains approximately constant 
when the position of the resonance moves through the protostar atmosphere. According to 
Eq. (j26j) ■ in the interval Rc ^ Rr ^ ^Vf^r allowed values of the oscillations parameters 
are 10^ eV'^ > Am^ cos26' > 10^ eV"^, that for small mixing leads to 30 KeV > m. > 1 KeV . 
From Eq. ()33|) it can be seen that transitions will be adiabatic for tan^ 26 > 7 x 10^^^ at 
Rr = Rc and tan^ 2^ > 3 x lO^^ at Rr = R- . 

In the second model, both the baryon de nsity and the temperature profiles follows simple 



potential laws and are related according to 



rj-] \ 3 



p = Pc ( ^ ) , (53) 



From this expression we see that = 3h^^, which in turn implies that hy^ = (see Eq. 
()3H) ). As a consequence, the electron fraction and also the ratio fie/T remain constants in 
the protostar atmosphere. Employing the same set of core parameters than in the preceding 
model we have = 0.08 and Pe/^ = 2.5. In what follows, for the density profile we use 
P = Pc (v^)^! which gives = — 4/r , and the border of the i/^ .^-neutrinoshere is located at 
p{Ru^,r) — 10^^ 9 cm~^, that corresponds to Rp^., ~ 3 i?c and T{Rp^^) ^ 9 MeV. The ratio 
of the diffusive term to the isotropic term in Eq. pn|) can be expressed as 10~^ (r/Rc)^'^^^ 
and the approximation (jl^ for A/C //C is also valid in this model. Taking these results into 
account and, from Eqs. ()29|1 and the following simple expression can be derived for the 
required magnetic field: 

5 = 3x10^^(^1 + 0.09^^ ^G, (54) 

with T = Tc{Rc/r)^^^. Now the intensity of the magnetic field decreases monotonically from 
3 X 10^^ G a.t Rr = Rc to 5 X 10^^ G at i?^ = Ruf^^- The allowed range for the sterile neutrino 
mass is similar to the one derived in the case of the Eddington model, while the requirement 
that the resonant conversion be adiabatic imposes the condition tan^ 2^^ > 3 x 10~^^ (Rr/Rc)^ 
on the mixing angle. 

A comment is in order. Our results have been derived on the basis of the weak field 
limit for the effective neutrino potentials in a magnetized medium. As it has been already 
mentioned in Section 2, this approximation is valid whenever B <C Pe/2e. In one of the 
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examples (the Eddington model) this constraint is verified in the interior but not in the 
external region of the protostar. A more careful treatment would require to incorporate the 
effect of strong magnetic fields along the lines examined in Ref. . 

V. FINAL REMARKS 

We have reexamined the anisotropic momentum emission by a nascent neutron star driven 
by the resonant neutrino conversion in presence of strong magnetic fields. The protostar 
atmosphere was modeled by a nonrelativistic nucleon gas and a degenerated electron gas. 
We elaborate the concept of a resonance spherical shell for neutrinos with a given momen- 
tum k moving in different directions relative to B. This concept provides a more precise 
description of the problem than a deformed spherical surface for neutrinos with an average 
radial momentum. An important ingredient for the existence of the effect is that the ingoing 
neutrinos become thermalized. For active-active neutrino oscillations {ue < — > ^h^t) the ther- 
malization proceeds by any of the following reasons: either because the ingoing z/^ moving 
toward the interior of the protoneutron star enter their own neutrinosphere, or because they 
cross again the resonance region and are reconverted into z/g, which are now within their 
neutrinosphere. For sterile neutrinos solely the second process is effective. In any case only 
the outgoing neutrinos contribute to the total momentum flux . The difference in the areas 
of the semispherical emission surface for neutrinos moving in opposite directions with the 
same momentum magnitude originates a non null fractional asymmetry in the total momen- 
tum. We apply our approach to the study of the generation of a natal pulsar kick through 
the vi < — > Us oscillations. The kick was calculated in a two simple model for the stellar 
atmosphere and in both cases a magnetic fleld of the order of 10^'' G is required to produce 
the observed velocities. The ranges of the mass and the mixing angle are compatible with 
a massive sterile neutrino as a warm dark matter candidate 122 ] . Fields of this magnitude 
are possible in protoneutron stars. For example, magnetic flelds of the order of 10^^ G have 
been considered in the literature for the central regions j^l, with a dependence in the baryon 



density parametrized as 



24| 



BpNS = Bs + Be 



(55) 
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with /3 = 10-5 and 7 = 3. In Eq. §^ Ps = W^^g cm-^, Bs = 10^^ G, and = 10^^ G are 
the strengths of the magnetic field at the surface and the core, respectively. For comparison, 
in Fig. 2 and 3 we have also shown the magnetic field as given by Eq. (|55|). Taking this 
curve as an upper value of the magnetic field in the protostar, we see that the resonant 
neutrino transformation could be an effective mechanism for the pulsar kick in most of the 
neutrinosphere. 
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Figure captions 

FIG. 1: Characteristic profiles for the Eddington model. The different functions are 
normalized to the corresponding value at the core: Rc — 10 km, Tc — 40 MeV, Pc — 
lO^'^g cm^, Mc = 1 Mq, and Fee = 0.08. 

FIG. 2: Magnetic field required to produce the kick, B, for the Eddington model and 
the expected maximum magnetic field, Bppfs^ functions of the resonance radius and 
normalized to the field required at the core, Be — 8 x 10^^ G. 

FIG. 3: Magnetic field required to produce the kick, B, for the potential model and 
the expected maximum magnetic field, Bpj^s, ^ functions of the resonance radius and 
normalized to the field required at the core. Be — 3 x 10^^ G. 
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